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Abstract 

We propose the estimates of the discrete Green function for the stream- 
line diffusion finite element method (SDFEM) on Shishkin meshes. 

1 Problem 

We consider the singularly perturbed boundary value problem 

(1.1a) Lu-~ -~eAu + b-\7u + u = f in O = (0, l) 2 , 

(1.1b) u = on dfi, 

where e <C 1 is a small positive parameter and b = {pi, &2) T > (0, 0) T is constant. 
It is also assumed that / is sufficiently smooth. 

2 The SDFEM on Shishkin meshes 
2.1 Shishkin meshes 

Let N > 4 be a positive even integer. We use a piecewise uniform mesh — a 
so-called Shishkin mesh — with N mesh intervals in both x— and y— direction 
which condenses in the layer regions. For this purpose we define the two mesh 
transition parameters 

1 „ e , „1 , , f 1 „ e 



\ x := min < -, 2— In N > and \ y := min < -, 2— In N 
{ 2 Pi J L 2 p 2 

Assumption 1. We assume in our analysis that e < N^ 1 , as is generally 
the case in practice. Furthermore we assume that \ x = 2e/3j~ 1 In TV and X y = 
2e/3^" 1 lnA r as otherwise N~ x is exponentially small compared with e. 
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Figure 1: Dissection of fi and triangulation Q N . 

The domain 51 is dissected into four parts as ft = Sl s U ft x U tt y U Sl^see 
FIG.Q}, where 

Q s := [0, 1 - X x ] x [0, 1 - A„] , n a := [1 - X x , 1] X [0, 1 - A. y ] , 
n y := [0, 1 - X x ] x [1 - X y , 1] , := [1 - X x , 1] X [1 - Ay, 1] . 

We introduce the set of mesh points {(xi, yj) £ SI : i, j = 0, ■ ■ • , iV} defined 

by 



and 



X j 



2i(l - A^/iV, for i = 0, • • • , N/2, 

1 - 2(iV - i)X x /N, fori = N/2 + 1, ••• , JV 

2j(l - A. y )/7V, for j = 0, • • • , iV/2, 

1 - 2(iV - j)V#, fOT j = N/2 + 1, ,N. 



By drawing lines through these mesh points parallel to the x-axis and y-axis the 
domain ft is partitioned into rectangles. This triangulation is denoted by ^^(see 
FIG. [I}. If D is a mesh subdomain of fl, we write D N for the triangulation of 
D. The mesh sizes h XtT = xi — Xi—\ and h yjT = yj — yj—i satisfy 



K~^, fovi = N/2 + l,-.>,N 
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and 

1 - A 



hy T — 



Hy:= ^]i> Sorj = l, — ,N/2, 
V=]^ for j = N/2+1,. ■■ ,N. 



The mesh sizes h X T and h l/ T satisfy 

AT 1 < H x , H y < 2N' 1 and deN^ 1 InN < h x ,h y < CaeA^ 1 InN, 

where C\ and Gi are positive constants and independent of e and of the mesh 
parameter N . The above properties are essential when inverse inequalities are 
applied in our later analysis. 

For the mesh elements we shall use two notations: r^- = .t.;] x [yj-i, yj] 

for a specific element, and r for a generic mesh rectangle. 

2.2 The streamline diffusion finite element method 

Let V := Hq(Q). On the above Shishkin mesh we define a finite element space 
V N := {/ e C(fi) : v N \ an = and v N \ T is bilinear, Vr e 
In this case, the SDFEM reads as 



(2.1) 



Find U e V N such that for all v N e V N 

e(VC7, Vw w ) + (b ■ VU + U, v N + 8b ■ Vv N ) = (/, v N + Sb ■ Vv N ). 



where 6 = 5(x) is a user-chosen parameter (see [3J). 
We set 



6 := + &2, /3 := f£j /&, ,7 : = f ^ 2 j /ft and « c := C T V V 

for any vector £ of unit length. By an easy calculation one shows that 
(Vw, Vu) = (W/J, + u,,). 

We rewrite (|2 . 1[) as 

e(U Pt v$) + e(U v , v%) + (bUp + U,v N + 6bv$) = (/, v N + 5bv$) 
and, following usual practice, we set 

0, otherwise. 



S(x) 



For technical reasons in the later analysis, we increase the crosswind diffusion(sec 
j3]) by replacing e(U v ,v^) by e(U v , v^) where 



max 



(e,N- 3 / 2 ) 



3 



and 

„, , _ J £j x £ fl s , 

We now state our streamline diffusion method with artificial crosswind: 
Find U £V N such that for all v N £ V N 

,N\ _ (f „,JV _i_ Xh,,N\ 



B(U,v N ) = (f,v» +8bv»), 



(2.2) 
with 

(2.3) B(U,v N ) := (e + b 2 6)(Uf,,v$) +e{U v ,tf) - 6(1 - 5)(U,vf) + (U,v N ). 

3 The discrete Green function 

Let x* be a mesh node in O. The discrete Green's function G £ V N associated 
with a;* is defined by 

B(v N ,G) = v N (x*),\/v N £V N . 
The weighted function w: 

uj{x) := ,g ff .g 

V a l3 J \ a v / V ^ 

where 

2 

g( r ) = ^— — r for r e ( _0O > °°)- 

1 + e r 

and cr^ = fcA^ 1 In AT and cr„ = ke 1/2 In AT. 
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= (6 + ^)11^/2^112+^-1/2^112 + J ||(w _ 1) 1/ 2g||2+||w _ 1/2g||2 



and 

(3.1) \\\G\\\l=B(u J - 1 G,G)-(s + b 2 6)((uj- l ) p G,G p ) 

-idu-^G^-bSiu-^Gp). 

Thus, we obtain 

G) = B(E, G) + Bduj^G) 1 , G) 
= B(E,G) + {uj- 1 G){x*) 

where E := uj- 1 G - (lu^G) 1 . 

Lemma 1. If erg = kN^ 1 In AT and a n = fclnA^e 1 / 2 , then for k > 1 sufficiently 
large and independent of N and e, we have 

B(u-iG,G)>±\\\G\\\l. 
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Proof. From (|3.ip ,wc estimate the following terms. 

(e + 8) KO^G, G P )\ < C(e + 5) l / 2 a~ p 1/2 ■ \\(u>-% /2 G\\ ■ (e + 5fl 2 \\^G j3 \\ 

V 

and 



<C{e + 5)^a- 1/2 \\\G\\\l 



e\{{LO- x ) v G,G v )\ < Ce^a- 1 ■ \\u>-V 2 G\\ ■ i^W^GJ 
<Ce 1 / 2 a-- 1 |||G|||2, 

For b5(uj~ 1 G, Gp), we make use of integration by parts. 

From the definition of ap and a v and e < N^ 1 , we take k sufficiently large 
and we are done. □ 

Lemma 2. If ap = kN~ IniV, with k > sufficiently large and independent 
of N and e. Then for each mesh point x* G f2 \ Q X y, we have 

K^G^x*)] < ±\\\G\\\l + CNlnN. 
where C is independent of N , e and x* . 

Proof. First let x* G £l s - Let r* be the unique triangle that has x* as its 
north-east corner. Then 



\(uj- 1 G)(x*)\ < CJV||G|| T . 



(u-% 1/2 ■\\(u-% 1/2 G\\ T , 



< CN max 

T* 

Calculating (oj^ 1 )p 1 (x) explicitly, we see that 

(uj- 1 )- 1 ^) < Cap = CkN' 1 InN G t* 

Thus 

\(u>-iG)(x*)\<CN]nN + ±\\\G\\\l 

by means of the arithmetic-geometric mean inequality. 

Next, let x* G ^.(The case x* G Vt y is similar.) Write x* = (xi,yj). Then 

lor 1 ^*)! = |G(x*)| 
1 

G x (t,y 3 )dt 



r 1 rVj+i 

<CH^ / \G x (t,y)\dydt 

Jxi Jy } 

^CNielnN-N- 1 ) 1 / 2 ^^ 
< CN 1 ' 2 ln 1/2 7V|||G||| 
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where G x {t,y 3 ) = Gk - ° k ~^ for (t, Vj ) £ T kj . 

Analysis: for the relation of boundary integral and domain integral, we analyze 



r 1 rvs+i 

Ixi Jyj 



N 



\G x (t,y)\dydt = h x 

k=i+l J Vi 



where f{y 3 ) = 



For Ai := Pf 1 ffa)*^ + f(y j+ i)*& 



_ Gkj+i—Gk-ij+i 



dy, we have 



Ai 



\f(yj)\ + \f(Vj+i)\ tt -1 r\fi \ \ft \i\tt 



if f(Vj)f(Vj+i) >0 



2 \f(Vi+l)-f(Vi)\ 
For Vw £ C 2 (t), we have 



> ^max{|/( % ), |/(y j+ i)|}ff„ if < 



□ 



Similarly, we have 



M <c(|« x | + K|) 

< C(\V XX \ + \l! X y\ + \Vyy\) 

Lemma 3. Let t e fl N . Then 

\\L0 1 / 2 D a E\\ Qs KCk-^N^WlGWl 
\\^ 2 D a E\\ QN \n B < Ck-'-e-^hT 1 N\\\G\\\ U 
where H T = in&x{h XtT ,h yjT } and \a\ = 1, D a G are Gp and G n . 
Proof. Assume p £ [1, oo] and g £ C 3 (r). Then (see [5J Theorem 4]) 

\\{g - 9 I )x\\lp{t) < G [hi !T \\g X xx\\LP(T) + h x>T hy jr \\g xxy \\ L p^ + h y>T \\g xyy \\ LP ^) 
+ Ch x ,T\\gxx\\Lp(T), 

11(5 _ 9 )v\\l>>(t) < G {h X T \\g xxy \\ L p( T ) + h x ^ T hy }T \\g xyv \\ LP ( t) + ^ !T ||.9 W !/||lp(t)) 
+ Ch y,rhyy\\Lp( T )- 

or (see [TJ Comment 2.15]) 

II (ff - 9 7 )x||L 2 (r) < C^,rll5a;ra||i 2 (T) + C^je.t ||fe||i 2 (r)) 
11(5 - 9 )v\\l 2 (t) < C^,tIIS^2/IU 2 (t) + Chy }T \\g yy || L2( T ) . 



G 



In the following analysis, D a v denotes the directional derivative of v along j3 or 
rj for different orders. The following analysis makes use of the former estimates 
(The latter will make the analysis more shorter). 
For t € Fl N , we have 

\\u 1/2 Ep\\ T < Cmaxu^iWExWr + \\E y \\ r ) 

T 

< Cmaxcj 1 / 2 -^ 2 T \\{u)- l G) xxx \\ r + h XtT \\{oj- l G) xx \\ T + h x , T H r \\(u)- l G) xxy \\ T 

T 

+ flr/h,,r||(w~ 1 G) ara ||r + ^.rlKW^ravllr + Kr II ("'^yy \\r} 

< Chl T {\\("~%™G\\T + \\(u-% x G x \\ T ) + Ch x , r {\\{uj-% X G\\ T + || (w- 1 )^,!,) 
+ ChlJW^yyyGWr + || (w^^G, || T ) + Ch y , T {\\ (aT 1 ) yy G\\ T + || (w^Gy || T ) 
+ C/ix^^dKa;- 1 )^^!^ + IKw-^xx^Hr + IKaT^GxIlT + ||(w _1 )xGxvI|t) 

+ flrftv,r(||(w~ 1 )«v»G ! ||T + IKW^^Ht + HOO^Hr + K^yGM 
3 2 

<CH*Y, \\D a (u- 1 )D"<G\\ T + CH T J2 \\D a (u- 1 )D"<G\\ T 

k = 2 |a| + | 7 |=3 k=1 |a| + |-y|=2 

o|>fc \a\>k 

where we have used the following analysis for J2\ a \=i M=2 D a (ui~ 1 )D~'G : 

\\^-%G x J T <C\\ E ^"(Ol-GxsHr 
|o|=l 

<Cmax V l^"^- 1 )! • \\G xy \\ r 

T Z — J 

\a\=X 

<Ch-*max £ l^^" 1 )! • ||Gy T 

|a|=l 

i«i=i 

The same analysis can be applied to ||a; 1 ' /2 £' n || T . 
For r S f2 s , we have 

\\lo 1/2 E p \\ t < Ck~ 5 / 2 N~ 2 [(a~ 5/2 + a/' 2 *- 1 + a/ /2 a- 2 )max(uj-% /2 + tr^maxu 
+ Ck-^ 2 H 2 [(a- 3/2 + af^a-^maxiu- 1 ) 1 / 2 + a^max^ 1 / 2 ] ■ N\\G\\ T 
+ Ck-V 2 H T [(<j- 3/2 + (j^^-Vax^- 1 )^ + tT-'maxw- 1 / 2 ] ||G|| T 

|a|=l | 7 |=1 

where we have used the estimates of w" 1 , standard inverse estimates and Holder 
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inequalities. Similarly, we have 

||« 1 /%|| n .<C* : - 1 /2jv 1 /2|||G||| (i; . 
For r G fl N \Of , we have 

||^ 1/2 ^||r < Cfc- 5 / 2 !^- 5 / 2 + a^V- 1 + ^ 1/2 a- 2 )max(^ 1 )y 2 + o^ma**" 1 /*] 
+ Cfc- 2 £- 1 / 2 iJ T 2 [a- 2 +a- 1 ( T- 1 +a- 2 ]max W - 1 / 2 .e 1 /2 £ || D 7 G || t 

ItI=i 

+ Ck- 3 ^H T [(a~ 3/2 + trJ^tT-^maxiu)- 1 ) 1 / 2 + u^maxw- 1 / 2 ] ||G|| r 
+ Cfc- 1 e- 1 / 2 i? r [^ 1 +(7- 1 ]max W - 1 / 2 . £ 1 / 2 ]T ||ZTG|| T 

l7l=l 

< C*- 1 ^- 1 / 2 lix -1 JV|||£3f|H„ 
Similarly, we have 

ll« V Xlln\n. ^ Cfc-^-^ln-^IHGIII^. 

□ 

Lemma 4. Lei r G Of . Let E ^ w _1 G - (w _1 G) 7 w/iere (o^G) 7 denote tte 
bilinear function that interpolates to lo~ 1 G at the vertices of t. Then 

||w 1/2 ^||n. ^Ck^N-^WGW^ 
\\^' 2 E\\ n \ a ,<Ck-^' 2 \\\G\\\ u 

where H T = max{/i IjT , hy,r\- 

Proof. We make use of the following standard interpolation error bounds 

||u - u'Wlp^) < hl T \\u xx \\ LP{T) + hl T \\u yy \\ LP(T) 

where p G [1, oo] and u G C(r) n W 2 ' p (t). 
Then, we have 

||£|| r < hl r \\{u- l G) xx \\ T + hl T \\{"- l G) V y\\r 

<CH 2 J2 p a (^ 1 )-G|| T + OT 2 ||(|( w - 1 ),| + |(c- 1 ),|).(|G /J | + |G„|)|| T 

\a\=2 

< CH? {Ku-^pGpWr + Kw-^pGJr + IKa;" 1 )^^ + \\iu-\GJr} 
+ CH 2 \\ J2 iD^co-^l-GWr. 

From the above inequality, we have 

\\u> 1/2 E\\n B < Ck-'N-^IWGW^ 
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and 



lk 1/2 £||n\n a < CArV 1 / 2 AT 1 |||G||| W . 
Following the techniques of (see [SI Lemma 4.4]), we have 

r r( a 



{0J 1/2 E)r,ds 



where x S f2 \ tt s , T(x) € T satisfies (a; — T(a;)) • (3 = and the following 
condition: 

For Vy G T, (x - y) ■ @ = 0, 

|a; — = niin|cc — y|. 

y 

From the above representation of oj x ^ 2 E, we have 



\\" 1/2 E\\l 



l-\y 

A 



i-a x Ji( a! „,r( a! „)) 



r(x) 



(w 1/2 £:)^ds 



dft 



(w 1/2 J B) I? ds 



1 2 






d«4 




[/ 




io Jl( Xtd .V( Xld )) 





r r( a 



an 



<cxl {11(^)^11^ + 11^/2^11^1 
^^in^l^ii^i^ + ii^/^ii^} 

< Ck~ 2 e 2 In 2 7V{iV 3/2 In" 2 N ■ e^N' 2 + e" 1 In" 2 N}\ \\G\ || 2 
<C£r 2 e|||G||| 2 



where Xq = j^X x and 

• € {(1 - X x ,y) : X < y < 1 - X y }; 

• ^zd G {(1 - K, v) ■ < y < A }; 

• x u € {(x, 1 — Aj,) : 1 — X x < x < 1}. 



□ 



Lemma 5. 1/(7,3 = kN 1 \nN and a v = ke 1 ! 2 In N , where k > 1 sufficiently 
large and independent of N and e. Then 

B{{^Gy-u- l G,G)<^\\\G\\\l. 
Proof. Cauchy-Schwarzs inequality gives 

\B(E,G)\ < (s + b'sy^Wu^EpW ■ (e + b 2 S) 1 / 2 \\cu^ 2 G p \\ + e^Wu^Ej ■ e^W^GJ 
+ C\\uj^ 2 E\\ ■ \\u-V 2 G p \\ + \\uj 1/2 E\\ • ||^- 1/2 G|| 

From Lemma [3] and Lemma 21 we are done. □ 
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Theorem 3.1. Assume that erg = kN^ 1 hiN and a v = ke 1 / 2 InN, where k > 
is sufficiently large and independent of e and N . Letx* £ Cl\Cl xy . Then for each 
nonnegative integer v, there exists a positive constant C = C{v) and K = K{v) 
such that 

\\G\\w^°°{a. B \n' ) < CN~ V , 
e\G\ W "i-.°o((n x un v )\n' ) + \\G\\L°°((a*uQ y )\ci' ) < CN' V 

and 

e\G\ w i,oo(a my \a^ + UGH^fo^n;,) < Ce' 1/2 N" V . 
Proof. On Q s , we apply an inverse estimate. 

On O \ tt s the application of an inverse estimate does not yield a satisfactory 
result, so we use a different technique. 

Let x E Q x \ tt' Q be arbitrary. Starting from x we choose a polygonal curve 
r C (O \ O x j/) \ t na t joints x with some point on outflow boundaries. If 
(x — x*) ■ r) < 0, we can choose T as a line parallel to f3. If (x — x*) ■ r\ > 0, the 
situation is a little complicated. We can choose T as follows: 
In fl x \ SIq, we choose the direction of T along r\ or the negative direction of 
x-axis so that T n £l xy = In (f2 s U n s ) \ Oq, we choose the direction of T 
along r\ or the positive direction of y— axis. 

Let T N be the set of mesh rectangle r in (f2 \ Vt xy ) \ that Y intersects. Note 
that the length of the segment of Y that lies in each r is at most CeN^ 1 In AT if 
r 6 fl x or r £ Q y . 

Then, by the fundamental theorem of calculus and inverse estimates in dif- 
ferent domain, we can obtain the results. □ 
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